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ABSTRACT 

A new catalogue of extrasolar planets is prese nted by re-analysing a selection of pub- 
lished radial velocity data sets using exofit ( Balan fc Lahav||2009| . All objects are 
treated on an equal footing within a Bayesian framework, to give orbital parameters 
for 94 exoplanetary systems. Model selection (between one- and two-planet solutions) 
is then performed, using both a visual flagging method and a standard chi-square 
analysis, with agreement between the two methods for 99% of the systems considered. 
The catalogue is to be made available online, and this 'proof of concept' study may 
be maintained and extended in the future to incorporate all systems with publicly 
available radial velocity data, as well as transit and microlensing data. 
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1 INTRODUCTION 

Since the discovery of the first extrasolar planet in 
1995 ( Mayor fc Queloz||1995[ ), the research on extrasolar 
planets has undergone exponential expansion. A wide range 
of search methods have been developed during this period, 
resulting in the discovery of more than 700 planets to date, 
the majority of which have been from the radial velocity 
method. Traditional data reduction methods use a peri- 
odogram ( |Lomb|[T976l |Scargle||i"982| ) to fix the orbit al pe- 
riod and then the Levenberg-Marquardt minimisation ( Lev- 
jenberg 1944 Marquardt 1963) to fit the other orbital pa- 
rameters. A catalogue of exoplanets has already been pub- 



lished by Butler et al. ( 2006 1 using this method to extract 
the orbital parameters of exoplanets. Recently, Bayesian 
MCMC me t hods have been introdu ced by |Gregory| ( |2005[ ); 



Ford ( 2005 I; Ford & Gregory ( 2007 1 as a replacement for the 



traditional data reduction pipeline, exofit (Balan & Lahav 
2009) is a freely available tool for estimating orbital param- 
eters of extrasolar planets from radial velocity data using a 
Bayesian framework. Here are analysed 94 previously pub- 
lished data sets using exofit, forming a new, uniformly de- 
rived catalogue of exoplanets from a Bayesian perspective. 

Statistical properties of the distribution of orbital pa- 
rameters are critical for explaining the planetary formation 
process. It has been argued that there is now a statistically 
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significant number of known companions to make inferences 
about the correlation between orbital elements. Early dis- 
cussions on this subject can be found in a series of articles 
on the statistical properties of exoplanets by |Udry et al.] 



(20031; Santos et al. ( 2003 1 ; EggcnbcrgeretaL (2004) and 



Halbwachs et al. (20051. The statistical discussion in this 



article is informed by the comparison to the published cat- 
alogues at http://www.exoplanet.eu ( Schneider et al.||2011 l 
and http://exoplanets.org (Wright et al. 20111. 

The rest of this article is structured as follows: in Sec- 
tions [2] and [3] the Bayesian framework and exofit software 
package are introduced. The data analysis pipeline is de- 
scribed in Section|4] model selection is discussed in Section[5] 
and the catalogue is presented in Section [6] and in the ta- 
bles at the end of the article. The statistical properties of 
the distributions of various orbital parameters are discussed 
in Section [7] and the results are summarised in Section [8] 



2 BAYESIAN FRAMEWORK 

The Bayesian framework provides a transparent way of mak- 
ing probabilistic inferences from data. It is based on Bayes' 
theorem, which states that for a given model H with a set 
of parameters © and data D, the posterior probability dis- 
tribution of parameters Pr(©|D,_ff) is proportional to the 
prior probability distribution Pr(0jff) times the likelihood 
of data, Pr(D|0, H). Using standard mathematical notation 
one can write, 



Pr(0|D,.ff) 



Pr(D|0,ff) Pr(0|H) 
Pr(D|#) ' 



(1) 
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The denominator of the right hand side of the above 
equation is called the Bayesian Evidence. Since it is the es- 
timation of parameters that is of interest here, this term can 
be considered as a normalising constant and Equation [l] can 
be written as, 



Pr(0|D,_ff) oc Pr(D|0,H)Pr(0|fl"). 



(2) 



The key step in the Bayesian approach is to obtain the 
posterior distribution of parameters accurately. The infer- 
ences are then derived from the posterior distribution. The 
Markov Chain Monte Carlo (MCMC) method is a widely 
employed technique for simulating the posterior distribu- 
tion (the left hand side of Equation [2|. The basic steps in 
Bayesian parameter estimation can be summarised as fol- 
lows: 

(i) model the observed data, i.e. construct the likelihood 
function, 

(ii) choose the prior probability distributions of parame- 
ters, 

(iii) obtain the posterior probability distribution, 

(iv) make inferences based on the posterior probability 
distribution. 

exofit is a software package that estimates the orbital 
parameters of extrasolar planets, following the steps outlined 
above. It should be noted that exofit does not perform any 
Bayesian model selection - for a discussion of the relation of 
this aspect of the Bayesian framework to this study, the 
reader is directed to Section [5] 



3 EXOFIT 



exofit (IBalan & Lahav 20091 is a publicly available tool 



for extracting orbital parameters of exoplanets from radial 
velocity measurements. It uses the MCMC method to simu- 
late the posterior probability distribution of the orbital pa- 
rameters. The likelihood of data Pr(Dj©, H) in Equation [2] 
connects the mathematical model to the observed data. The 
radial velocity model and the corresponding Gaussian like- 



lihood function are given in Balan & Lahav (20091, where 
the choice of likelihood is based on |Gregory (|2005 |. The 
prior probabilities are as used by Ford & Gregory (20071. 



exofit then generates samples from the posterior distribu- 
tions of the orbital parameters in the mathematical model, 
which can be analysed with the aid of any statistical soft- 
ware. Details of the algorithmic structure of the code, in- 
cluding methods of controlling chain mixing and assessing 



convergence, can be found in Balan & Lahav ( 2009 1. For fur- 



ther information, the reader is directed to the exofit user's 
guide. 



4 DATA ANALYSIS 

As of 2009 August 21, when the data sets were ex- 



tracted from the literature (Butler et al. (2006), and ref- 
erences therein), there were 295 planetary systems de- 
tected using the radial velocity method according to 
http://www.exoplanet.eu, with 346 individual planets in to- 
tal. Some radial velocity data sets from the literature have 
fewer than ten entries and as such are not appropriate for 



use with Bayesian inference methods, as a radial velocity 
data set will need to include at least half an orbital period 
of a potential planetary companion. Hence data sets with 
not enough measurements to give accurate orbital solutions 
were not included. Also, the radial velocity data of any sys- 
tems with more than 2 confirmed planets were ignored since 
at present, exofit can only search for either one or two 
planets. 

Many more and different radial velocity data sets and 
stellar mass estimates are now available (though some not 
publicly), but for the sake of uniformity the original radial 
velocity data were used (i.e. those publicly available, frozen 
as of 2009 August 21 when the original data were collected). 
At a later stage the results can be improved by updating 
the original data sets to those which are now available, as 
well as incorporating the data from the many hundreds of 
additional planets that have been detected since the start of 
this study. 

To enable accurate calculation of the derivable orbital 
parameters, the masses as well as the radial velocities of the 
associated stars were needed. These values were all taken 
from the published literature at http://www.exoplanet.eu, 
frozen as of 2011 March 01. The input for exofit is in the 
form of a simple text file with radial velocity, uncertainty, 
and the time of observation (in Julian Date format), where 
the radial velocity values must be in ms" 1 . The Julian Date 
of the observation is offset to zero within exofit. 

The publicly available statistical data analysis pack- 
age, R, from the R Project for Statistical Computing 
(http://www.r-project.org) , was used to analyse the output 
of exofit. The output from R includes the mean, median 
and standard deviation of the orbital parameters extracted 
from the posterior distribution samples produced by exofit. 
The modal values are also produced, but will only have sig- 
nificance in the event of the posterior having more than one 
peak. Posterior distribution plots can also be produced with 
R, and the marginal distributions of each parameter can be 
found by plotting a histogram of the samples from the pos- 
terior. The full posterior distribution is helpful in analysing 
correlations between various parameters. Even though pa- 
rameter degeneracy is present in the orbital solutions, highly 
degenerate solutions are less common. 

The calculation time required for exofit depends on 
computational resources available to the user. It scales lin- 
early with the number of radial velocity entries input to 
the code, and also depends significantly on the ease with 
which exofit can converge the data. If exofit is pre- 
sented with data from a non-converging posterior distribu- 
tion, it will take much longer than a larger data set with 
convergent orbital solutions. In technical terms, the mean 
average calculation time using exofit on 26 radial veloc- 
ity data sets ranging from 10 to 50 data entries for a 1- 
planet search was 44 seconds per data entry. For a 2-planet 
search using 30 data sets with between 11 and 256 entries, 
the calculation time increased to 3 minutes and 40 sec- 
onds per data entry. These times were achieved on a 2.80 
GHz dual core linux system. Multiple runs were performed 
in order to confirm the orbital solutions for each system 
- these analyses were carried out using the UCL Legion 
High Performance Computing Facility, details of which can 
be found at http://www.ucl.ac.uk/isd/common/research- 
computing/ services /legion-upgrade. 
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5 MODEL SELECTION 

One of the most challenging aspects of the statistical in- 
ference procedure is the model selection problem. For the 
analysis of the radial velocity data the question of model se- 
lection refers to the selection of the correct number of plan- 



ets to fit the observed data. Ford & Gregory ( 2007 1; Gregory 
( 2007 1 employed thermodynamic integration for calculating 



the Evidence and selecting the optimal number of planets 

( j2011b|a l 



that fit the data. On the other hand, Feroz et al 



approached the situation as an object detection problem. 

One of the most commonly employed model selection 
procedures makes use of the chi-square statistic. This is one 
of the most prominent methods for estimating the goodness 
of fit and it has been applied to many astronomical problems 



including the analysis of radial velocity data (sec e.g. Butler 
et al. | |2006| )). Bayesian inference also offers a straightfor- 
ward way of performing statistical model selection, based on 
Equation [T] and the Evidence. Even though this approach is 
conceptually simple, its implementation is in general com- 
putationally expensive, and exofit does not currently have 
the functionality to perform such Bayesian model selection. 

Hence, in the present analysis, we make use of the tra- 
ditional chi-square statistic as well as a visual flagging ap- 
proach, discussing the relationship between the two in Sec- 
tion]^] We limit ourselves to 1 or 2 planets, as per the current 
capabilities of the code, but this may of course be extended 
in later studies. The rationale behind the visual flagging is 
that one can identify the poor fits to the data by comparing 
the predicted radial velocity curves for the 1-planet and the 
2-planet solutions. The method involves assigning a 'visual 
quality flag' by eye to each system, where '1' signifies that 
the 1-planet fit is best, '2' means that the 2-planct fit is 
best, and '3' means that both 1- and 2-planet solutions pro- 
vide equally good (or equally bad) fits. The results of this 
classification are shown in Table |6j next to the number of 
planets currently confirmed to exist in that system, taken 
by comparing the values on both http://www.exoplanet.eu 
and http://exoplanets.org (as these catalogues do not agree 
with each other in some cases). 

Table[6]also shows the log likelihood ratio of the reduced 
chi-square value of the 1-planet fit to that of the 2-planet 
fit, where we define the log likelihood ratio, 



R 



(Xip 



2 \ 
X2p) 



(3) 



Hence a value of R > indicates that the 1-planet fit 
was best (had a smaller reduced chi-square value) , and R < 
indicates that the 2-planet model provided the best fit to 
the data. For all bar one system (HD8574), every '1' and 
'2' quality flag assigned to the fits by eye was in agreement 
with the calculated value of R, endorsing our method of 
assignation by visual inspection (see Figure [5|. For only a 
few systems there were not sufficient degrees of freedom to 
calculate a value for R (due to e.g. only having 11 datapoints 
for the 12 parameters), denoted by '-' in the table. 



6 CATALOGUE OF EXTRASOLAR PLANETS 

In this paper the catalogue of extrasolar planets generated 
using exofit is presented in Tables [4] and [5] These contain 
the best estimates of the orbital parameters for both 1- and 



2-planet fits for all systems analysed. The orbital parameters 
used to fit to the model were the systematic velocity offset 
of the data V, the orbital period of the planet T, the radial 
velocity semi-amplitude K, the orbital eccentricity e, the ar- 
gument of periastron w, and x> parameterising the fraction 
of the orbit at the start time of the data along which the 
planet has travelled from the point of periastron passage. 
The final parameter, s, is a measure of all extra signal in 
the data after the planetary fits have been accounted for, 
and hence a high value could indicate the presence of an 
additional planet, or noisy data due to stellar activity, or 
the combined noise from all sources. The reader is referred 
to |Balan fc Lahav ( 2009 1 for a more complete description of 



this parameter, which is not considered in any more detail 
in this study. It should also be noted that the orbital pa- 
rameter x is n °t related to the statistical measure \ 2 used 
for determining the log likelihood ratio in Equation [3] 

The direct output values from exofit shown in the 
tables are the medians of the parameter posterior distri- 
butions, and the associated 68.3% confidence regions. The 
other displayed derivable parameters of the systems (mass 
and semi-major axis) were calculated by transforming the 
orbital parameter posteriors using the standard relations, 



and, 



m p sim 



(2ttG) 1 /3 



(4) 



(5) 



assuming m p <C m». The final values for these deriv- 
able quantities were again taken to be the medians with 
68.3% confidence regions, and are also displayed in the ta- 
bles. 



6.1 Choice of priors 

The prior distributions and ranges used for the initial anal- 
ysis were as shown in Tables [I] and [2] The prior for the sys- 
tematic velocity was dependent on the system - the mean 
of the input RV data was calculated and used as the initial 
value, and the allowed range was 10 fcms -1 symmetrically 
about this. For some systems different sets of prior bound- 
aries were used in a second round of analysis - these stars 
and the prior boundaries applied are listed in Table [3] Sys- 
tems which did not return good fits using the normal prior 
boundaries were re-run with these 'tight' priors, where the 
period of the planet was constrained to be within a range 
given by, 



(6) 



where T pu b is the published value of the period and 
o~ P ub is the published error on the period, both taken from 
http://www.exoplanet.eu on 2011 August 01. Further con- 
straints may also be applied, for example, systems with near 
zero eccentricities require tight priors on the orbital param- 
eter x m order to avoid multimodal distributions (see Sec- 
tion]^}, whilst systems with eccentricities close to unity need 
tight priors on the orbital period T in order to achieve con- 
vergence of the MCMC chains. Examples of the output of 
exofit are shown in Figures [l] and [2] 
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Table 1. The assumed orbital parameter prior distributions and 
their boundaries for a 1-planet model. The min and max values 
for the systematic velocity parameter were the mean value of the 
raw radial velocities for that data file minus and plus 5000ms -1 
respectively. 



Parameter 



Prior 



Mathematical form 



Min Max 



Todays) 



Uniform 
Jeffreys 



Ki(ms 1 ) Mod. Jeffreys 



(K l +K l0 ) 



K 1 g + K 1 



0.2 15000 



0.0 2000 



<=1 
Xl 



Uniform 
Uniform 
Uniform 

Mod. Jeffreys 



1 



1 

2tT 
1 



Table 2. The assumed orbital parameter prior distributions and 
their boundaries for a 2-planet model. The boundaries for V were 
as detailed previously. 



Para. 



Prior 



Mathematical form 



Todays) 



Uniform 
Jeffreys 



Kx(ms - 1 ) Mod. Jeffreys 



e.\ Uniform 

uj\ Uniform 

Xi Uniform 

T% (days) Jeffreys 



K 2 (ms x ) Mod. Jeffreys 



(K 1 +K 1Q ) 



( Ki + K lr 
V K l 



1 

1 

2;r 
1 

JL 

^2 , 



(J 



0.2 



1 
1 

15000 



(K 2 + K 2() ^ 



£2 

X2 
s(ms~ 1 ) 



Uniform 
Uniform 
Uniform 

Mod. Jeffreys 



( K 20 + K 2max \ 

\ K20 J 



1 n 
1 

(«+° r 



1 

2tt 
1 



Table 3. Radial velocity data sets analysed with different period 
prior boundaries, for reasons explained in Section |6.2| The initial 
value is set to the published value of the period, the maximum 
value is the initial value plus twice the published error, and the 
minimum is the initial value minus twice the published error. This 
approach was generally necessary for those systems (e.g. WASP 
and XO data sets) where the number of datapoints available at 
the time of selecting the data was low, thus requiring tighter priors 
to adequately constrain the solution. 



System 



Initial period value Min period Max period 



cpsilon Eri 


2500 


1800 


3200 


gamma Ccp 


906 


899.84 


912.1 


GJ849 


1900 


1400 


2400 


GJS6 


15.7649 


15.76412 


15.76568 


HAT-P-9 


3.92289 


3.92281 


3.92297 


HD118203 


6.1335 


6.1323 


6.1347 


HD12661 


262.71 


262.54 


262.88 


HD128311 


924 


913.4 


934.6 


HD131664 


1950 


1868 


2032 


HD142 


350 


342.8 


357.2 


HD149143 


4.072 


4.058 


4.086 


HD162020 


8.42820 


8.428088 


8.428312 


HD168443 


58.1121 


58.111142 


58.113058 


HD169830 


225.6 


225.16 


226.04 


HD183263 


627 


624.8 


629.2 


HD187123 


3.096583 


3.09656732 


3.09659868 


HD189733 


2.2185757 


2.2185754 


2.2185760 


HD190360 


2920 


2862.2 


2977.8 


HD196885 


1330 


1300 


1360 


HD202206 


255.87 


255.75 


255.99 


HD20868 


380.85 


380.67 


381.03 


HD209458 


3.5247486 


3.52474784 


3.52474936 


HD217107 


7.12682 


7.1267318 


7.1268882 


HD219828 


3.833 


3.807 


3.859 


HD28185 


379 


375 


383 


HD330075 


3.38773 


3.38757 


3.38789 


HD33636 


2128 


2111.6 


2144.4 


HD38529 


2146 


2134.98 


2157.02 


HD46375 


3.02357 


3.02344 


3.0237 


HD47536 


430 


0.2 


860 


HD50499 


2460 


2384.2 


2535.8 


HD5319 


670 


636 


704 


HD68988 


6.2771 


6.27668 


6.27752 


HD73267 


1260 


1246 


1274 


HD74156 


2520 


2490 


2550 


HD75289 


3.50927 


3.509172 


3.509398 


HD80606 


111.4367 


111.4359 


111.4375 


HD86081 


2.1375 


2.1371 


2.1379 


HD89307 


2170 


2094 


2246 


HIP 75458 


511.10 


510.922 


511.278 


tau Boo 


3.31246 


3.312432 


3.312488 


TrBS-3 


1.31 


0.2 


2.62 


WASP-2 


2.152226 


2.152218 


2.152234 


WASP-3 


1.846834 


1.84683 


1.846838 


XO-1 


3.94153 


3.941476 


3.941584 


XO-2 


2.615838 


2.615822 


2.615854 


XO-4 


4.125083 


4.125075 


4.125091 



6.2 Ambiguous systems 

For some of the systems analysed there is a clear trend in 
the radial velocities indicating the possibility of a second 
planet, but the data are not informative enough to properly 
constrain the orbital parameters of such an object. Plotting 
the resulting radial velocity curve and judging by eye can 
help to assess and distinguish between the 1- and 2-planet 
fits and evaluate the validity of the orbital solution, though 
such poorly-constrained orbits will lead to large error bars 
on the estimates of the orbital parameters. 

This did not always even lead to a clear classification 



though, and many systems were then re-run with tight priors 
on the period, given in Table [3] Some of these ambiguities 
were caused by data which were poor, or less accurate due 
to age, or too noisy due to stellar jitter. Others were simply 
due to the correlation between u and x> or data not good 
enough to constrain these two parameters. This resulted in 
near-uniform posteriors for uj and Xi an d hence fits that 
match few of the measured datapoints as a result of being 
shifted in time. Estimates for masses and semi-major axes 
derived from these results are still valid however (provid- 
ing reliable estimates for T, K and e are obtained, which 
was generally the case), as these values have no dependence 
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Figure 1. The resulting marginal posterior distributions of the 
orbital parameters for a 1-planet fit to the BD-f7 63 data, and 
the corresponding radial velocity curve. 



Figure 2. The resulting marginal posterior distributions of the 
orbital parameters for a 2-planet fit to the HDf08874 data, and 
the corresponding radial velocity curve. 




(a) Posterior distributions of orbital parameters for BD-17 
63. 




53000 53500 54000 545 

Time (days) 

(b) Radial velocity plot for BD-17 63. 



on mean anomaly at epoch and the time evolution of the 
Keplerian orbit. 

The class 3 (both 1- and 2-planet fits equally good, or 
equally bad) systems, as introduced in Section [5] are those 
which were considered to be somewhat ambiguous even af- 
ter being analysed with tighter priors. This category was 
subdivided further - in some cases these are distinct radial 
velocity solutions which provide plausible fits for both 1- 
and 2-planet models, and are classified as '3a'. However, 
there are also systems where the 'second planet' fit just pro- 
duces small-amplitude variations on the 1-planet solution 
(see Figure [3] for an example), or where the 1- and 2-planet 
fits are identical but the 'second planet' posteriors are peaks 



14 16 18 20 22 1500 1600 1700 1800 1900 

(a) Posterior distributions of orbital parameters for 
HD108874. 




11500 12000 12500 13000 13500 

Time (days) 



(b) Radial velocity plot for HD108874. 

at very small values for T and K, and uniform for e, w and 
X (i.e. there is no single solution for a second planet from 
these data). From this an 'Occam's Razor' approach could 
be taken and the assumption made that the correct model 
for most of these '3b' class systems is in fact the single planet 
one. In a few cases though there may truly be a second planet 
present, and the data used are simply not good enough to 
change the likelihoods of the parameters from the initial 
'no-knowledge' (uniform prior) situation. So for all class 3 
systems, better (or at least more up-to-date) data and more 
complete analyses (such as using the log likelihood ratio in 
more detail to narrow down the classification) are required 
to accurately determine the correct orbital solution. 
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Figure 3. The resulting radial velocity curve fits using the de- 
rived orbital parameters for HD89307. Imposing a 2-planet model 
on data with only one planet can have the result shown here, 
where the 2-planet fit is the same as the 1-planet fit with a su- 
perimposed artificial small-amplitude periodic variation. 




11000 



12000 



11500 
Time (days) 

(a) 1-planet radial velocity plot of HD89307. 



^ o — 




11000 



12000 



11500 
Time (days) 

(b) 2-planet radial velocity plot of HD89307. 

7 COMPARISON OF ORBITAL PARAMETERS 

Figure [4] shows values for specified orbital elements from the 
literature against values yielded using exofit for each sys- 
tem. The mass, semi-major axis and period values all exhibit 
good correlations in general between the independently de- 
rived values and those in the published literature - this is 
unsurprising for the period as it has not been derived from 
other quantities. The eccentricity however shows a greater 
spread than expected - whereas our uniformly derived cata- 
logue is consistent in taking the median, the published values 
use, in general, many different statistical measures to extract 



the final parameter values from varying analysis techniques. 
This highlights the value of using a consistent technique to 
build up reliable databases of orbital parameters. 

As mass and semi-major axis values are themselves de- 
rived from period and eccentricity, any inaccuracies in algo- 
rithms used will propagate, and also present discrepancies 
in the values yielded using exofit and are likely to amplify 
outliers in these plots. These outliers will be investigated in 
the future in order to assess the validity of the solutions. 

There are some discrepancies in the global distribution 
of parameter values between this catalogue and the pub- 
lished literature, especially for the eccentricity parameter. 
This may be partly due to poor or out-dated data, and is 
almost certainly affected by the ubiquitous effects of certain 
parameter correlations (as explained in Section [81. These 
should be analysed in more detail in the future, and tech- 
niques developed to explore the parameter space more effi- 
ciently and minimize or eradicate such dependencies. 



8 DISCUSSION 

The primary objective of this article is to analyse radial ve- 
locity data sets uniformly, using a single platform for the 
data analysis. Butler et al. (20061 produced a catalogue 

(using pe- 
We 



of extrasolar planets using traditional methods 
riodograms and Levenberg-Marquardt minimisation), 
have analysed a selection of radial velocity data sets us- 
ing a Bayesian parameter estimation program for extrasolar 
planets. However, a model selection criteria is required for 
completion of the statistical inference process, and for this 
purpose, as described in Section[5] a chi-square statistic was 
employed as well as a visual flagging technique. Inconclusive 
results are obtained for a few data sets, but from those anal- 
ysed here it can be seen that both model selection methods 
perform well, agreeing in 99% of cases, as demonstrated in 
Figure [5] 

Investigating further, we find that the chi-square val- 
ues depend on the point estimates of the orbital parame- 
ters used to construct the predicted radial velocity curve. If 
the posterior distribution is unimodal such an approach will 
work flawlessly. However, posterior distributions of the or- 
bital parameters exhibit multi-modality on many occasions. 
For example the parameters uj and \ are extremely corre- 
lated and their posterior distributions are bimodal for many 
data sets (an example of this is shown in Figure[6|, especially 
for planets with e w 0. This problem has also been noted 



by Gregory (20071, who proposed re-parameterisation of the 



problem as a possible way of dealing with this situation. 

Many data sets contain planetary signals whose period 
is greater than the span of the observations, and so obtaining 
constraints on the orbital parameters of these objects is an 
extremely difficult task. There are several data sets where it 
was possible to obtain estimates for the orbital parameters 
for one of the planets, but then the second signal could not 
be constrained due to weak signal-to-noise. In most cases 
these signals appear to be a linear or quadratic trend in 
the radial velocity data. Therefore, it becomes extremely 
difficult to classify these objects as planets, and this is one 
of the reasons a visual flagging method was employed. One 
example of this is shown in Figure [6j the results of the 2- 
planct fit to the data of the system HD190228. The strongest 
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Figure 4. Orbital parameter values taken from http://www.exoplanet.eu, plotted against values yielded using EXOFIT. Plotted systems 
are only those where EXOFIT gave unambiguous (either class 1 or 2) results. 
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(c) Mass (in M Jup ) 
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Figure 5. The log likelihood ratio for each planet assigned to visual class 1 or 2 is shown in |( a) | [(b)] shows the same data on a smaller 
scale, around the threshold at R = 0. Class 2 systems (open circles) are all below the chi-square ambiguity threshold, and class 1 systems 
(filled circles) are all above, with the single exception of HD8574 (shown as a red triangle), class 1 but located just below the threshold 
with a value of R = —0.03. 
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(a) Log likelihood ratio with planet for all class 1 and 2 
systems. 



(b) Log likelihood ratio with planet for all class 1 and 2 
systems, in the range —20 < R < +20. 
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Figure 6. Examples of bimodal and ambiguous posterior densi- 
ties, obtained from HD49674 and HD190228. 



We have also found that the aliasing effects (see e.g. Dawson 
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(a) Posterior densities for a 1-planct fit to the HD49674 data, 
exhibiting some bimodality in the lj and \ orbital parameter 
values. 



& Fabrycky (2010 1) in observations can produce additional 



peaks in the posterior distributions, necessitating the use of 
the sharp prior on the period. 

In summary, a brief overview of the Bayesian theory 
has been given here, along with a description of the MCMC 
approach used in order to estimate the orbital parameters 
of extrasolar planets, more details of which can be found 
in | Balan fc Lahav| ( 2009 1. A new catalogue of extrasolar 



planets is presented from the re-analysis of published ra- 
dial velocity data sets, giving both 1- and 2-planet orbital 
solutions for 94 systems derived on a uniform basis. An at- 
tempt is made to distinguish between the solutions for each 
system by using both a visual categorisation method and a 
standard reduced chi-square technique, giving good agree- 
ment in 99% of the cases presented here. Improvements in 
this 'model selection' area of the analysis may be made by 
taking into account Bayesian Evidence, as seen in [Gregory] 
( [20071 ) and |Feroz etaLl ( |2011b|a| ); more rigorous approaches 
such as these are outside the scope of this 'proof of concept' 
study, but may be looked into in the future. Other further 
work will include updating this catalogue to incorporate the 
most up-to-date data, as well as extending exofit to be 
able to use transit and microlensing results, to search for an 
arbitrary number of planets, and to look into the possibility 
of accounting for interactions between planetary bodies. 
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(b) Posterior densities for a 2-planet fit to the HD190228 
data, showing ambiguity in the estimates for the ui, \, and 
e values for one of the planets. 



signal is picked up and well-constrained, as can be seen from 
the error bars in Table [5j and in addition the values for 
this planet match well those from the fit for a single planet. 
Thus we can be reasonably sure of the parameters of the 
first planet, but those of the second, shown as HD190228b 
in Table [5j are significantly less secure. 

Additionally, sharp prior boundaries were used on the 
orbital period for several data sets. In these cases we have 
found that either the planetary signal is very weak or the 
signal from a systematic trend from an additional companion 
in the radial velocities masks the weaker planetary signal. 
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Table 4. Table of the orbital parameters for a 1-planct fit, both directly output from EXOFIT and thence derived. The values of the 
parameters T, K, e and s (generated from exofit) are the medians of the parameter posterior distributions, with the associated 68.3% 
confidence regions. The other parameters were calculated using these values and stellar masses taken from the published literature. Note 
that some parameters are extremely well-constrained, hence the errors on the parameter estimates are so small as to appear to be zero 
to the two decimal places shown in this table. A full table in machine-readable format will be provided on the website, and the reader is 
directed there if such data arc required. 
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29+ 00 
u -^ a -o.oo 

,+0.01 



1.03 
0.07 
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planet m* (M soi ) T(days) K(ms 1 ) e s m p (Mj up ) a(AU) 



HD4203 


b 


1 


13 


434.21+^; 6 / 7 


74 24+ 88-38 
-18.24 


0. 


71 t°o 


.09 


5 R3+ 1 - 49 


2. 


.15 


+ 1 
— 


42 
.39 


1 


.17+.°, 


00 


HD4208 


b 





.88 


828.90+^°! 


19 01+°- 68 


0. 


05+0 


.03 


1 - 2b -0.85 


0. 


81 


+0 
-0 


.03 
.03 


1 


■ 65 -o 


.01 


HD43691 


b 


1 


.38 


36.93+.°,;°* 


123 


7 „+4.76 
78 -4.65 


0. 


"to 






2. 


.50 


+0 
-0 


.10 
.10 





- 24 lo 




HD43848 


b 





.93 


2390.59+^ 1 Q | 7 


888 


, 7 +447.92 
-286.53 


0. 


81+°- 


06 
.08 


5 - 58 -2-.90 


33.07+°;^ 


3 




11 

.08 


HD46375 


b 





93 


r> no+0-00 

a - u -o.oo 


•"■°'_0.80 


0. 


06+° 




.> n + ■ 59 
■ SiU -0.52 


0. 


.23 


+0 
-0 


.01 
.01 





.04+0 


.00 


HD47536 


b 





.94 


717.27+ J!;™ 


108 


71 +12.96 
-12.35 


0. 


!5l°o 


10 

.09 


6 65 + 9 ' 40 


4 


.53 


+0 
— 


57 
.56 


1 


-54+.0 1 


02 


HD49674 


b 


1 


.01 


4 qe+0.00 

*- ao -o.oo 


11 76+ 1 ' 09 


0. 


05t°0 




3 65 + 0- 78 
d0t, -0.69 


0. 


.10 


+0 
-0 


01 
.01 





.06±° 


.00 


HD50499 


b 


1 


28 


2460.46+^°° 


ob 77 +5.54 
ij ' 77 -4.56 


0. 


29+° 


9 1 
.19 


13.38+ 2 , 8 


1 


91 


+0 
-0 


.33 
.32 


3 




.05 


HD5319 


b 


1 


.60 


684.51^;^ 


36.23+ 1. 7 8 ° 9 


0. 


io+° 


20 


10 44+ 1 ' 73 
-1.37 


2. 


.13 


+0 
-0 


91 
.35 


1 




02 
.03 


HD63454 


b 





.80 


2.82+°°° 


63 32+ 1 '' 76 


0. 


02+.°, 


03 


5 - 8 °to.98 


0. 


38 


+0 
-0 


01 

.01 





.04+0 


.00 


HD68988 


1 


12 


6.28t°° ° 


183 


+ 13.27 
-13.39 


0. 


09+.°, 


or 

.06 


42.29l 7 ; 7 8 


1 


79 


+0 
-0 


12 
.13 





07l° 


.00 


HD73267 


b 





.89 


1259.62+°;" 


RA 9O + 0.44 
64 ' 28 -0.46 


0. 


26+° 


01 


°- 72 to. 5 48 


3. 


06 


+0 
— 


02 
.02 


2 


.20+.°, 


.01 


HD73526 


b 


1 


.01 


193.32+°-^ 


114 


7q + 22.33 
-18.22 





57+.°, 


09 




2. 


.70 


+0 
-0 


-10 
.35 





.66+0 


.00 


HD74156 


b 


1 


.24 


2519.62+ 2 2 °;f 3 


120 


„„ + 139.67 
88 -42.57 


0. 


.89+.°, 


.13 




4 


31 


+1 

-0 


34 
.68 


3 


.89+0 


.02 


HD75289 


b 


1 


.19 


3 - 51 -o°o°o 


53.90+;-;^ 


0. 


Oilo 


01 


o.73+i;° 4 


0. 


45 


+0 
-0 


.01 
.01 





™t° 


.00 


HD75898 


b 


1 


.30 


421.64+™ 


74. £ 


IQ +3.83 
-3.66 


0. 


05+°- 




9.56+ 2 ;f g 


3. 


29 


+0 
-0 


.17 
.16 


1 


.20+0 


.01 


HD76700 


b 


1 


13 


3 Q7+O.OO 

,ia7 -o.oo 


27.30+j'^ 


0. 


08+° 


06 
05 


2.79+i-f, 


0. 


23 


+0 
— 


01 
.01 







00 


HD80606 


b 





.96 


111 44+°-°° 
-0.00 


560 


6 „ + 187.15 
" -113.26 


0. 


95+° 


02 


i2.47+. 2 ;f 3 


3. 


84 


+0 
— 


.51 
.30 





•«ts- 


.00 


HD81040 


b 





.96 


1108.32+^;^ 


169 


71 +10.03 
-10.03 


0. 


42 +0. 


.14 


25 - 51 ±5.16 


7. 


62 


+0 
-0 


56 
.47 


2 




.01 


HD82943 


b 


1 


.13 


4030.24+™^- 5 2 2 3 


88.77+ 8 f 1 92 


0. 


49+° 


.33 


34 40+ 8,27 
° -6.74 


5. 


.58 


+ 28.20 
-4.69 


5 


16+ 4 




HD8574 


b 


1 


12 


227.45+° ; 8 7 2 9 


64.92+ 4 -f T 


0. 


28+° 


.05 


9 70+ 2,53 
a - -2.39 


2. 


02 


+0 
— 


.13 
.13 





- 76 ±o 


00 
.00 


HD86081 


b 


1 


.21 


214 to.°oo 


207 


49+0-86 
-0.84 


0. 


oi+° 


01 


76+0- 92 
u ' -0.56 


1 


49 


+0 
— 


01 

.01 





.03+0 


.00 


HD89307 





.99 


2174.23^f 3 


38.05+ 2 « f 


0. 


29+° 


.20 


r> +5.42 
J " 31 -2.68 


2. 


29 


+0 
— 


83 
.41 


3 


- 27 lo 


.05 


HD93083 


b 





.70 


143.99tJ; 47 


18. t 


■ 0+1-26 
lS -1.20 


0. 


i 2 io 


.06 


2.28+.0-75 


0. 


38 


+0 
-0 


.03 
.03 







00 


HR 810 


b 


1 


.11 


302.94+ 2 ; 27 


57 17+5-60 
j7 - 17 -5.43 


0. 


12 ^0 


08 


17.92+. 3 ;™ 


2 


00 


+0 
-0 


20 
.20 





91 +° 

al -o 


00 


kappa CrB 


b 


1 


80 


1218. 78+ 36 - 22 
— 28.41 


24.06+J' 89 


0. 


12+°- 

-0 


08 
.07 


4 73+0.86 
-0.75 


1. 


86 


+0 
-0 


.11 

.11 


2 


.72+0 
-0 


05 
. 04 


NGC 2423 3 


b 


2 


.40 


713 95+ 4 - 99 


133 


14 + 7-95 
-7.28 


0. 


l«+° 


06 
.06 


18.05+ 2 ; 72 


1 32 + 0-56 
iUJ ' ! -0.55 


2 


09+0 


01 

.01 


NGC 4349 127 


b 


3 


.90 


676.93+*'^ 


189 


+ 10.22 
^*°-9.44 


0. 


l^°o 


05 
05 


14 59+ 3 43 


19 89+ 1 ' 06 
a ' oa -1.02 


2 


,38+0 


01 
.01 


tau Boo 


b 


1 


34 


+10.38 
J ' 3 -10.09 


469 


„ 5 + 14.76 
ut> _14.87 


0. 


orlS- 


03 
04 


96.65+ 8 - 42 


4 


17 


+0 
— 


.13 
.13 







00 
.00 


TrES-3 


b 





.92 


1 34+ - 04 
-0.04 


352 


oe + 17.08 
-21.08 


0. 


06t° 


09 
.04 


10 70+ 24 41 
-9.07 


1 


80 


+0 
— 


09 
.12 





,02+0 


00 

.00 


WASP-2 


b 





.88 


9 . k+0.00 
j!10 -0.00 


159 


43+8.00 
-6.39 


0. 


25+° 


10 
11 


9 „ 9 + 5.14 


0. 


.90 


+0 
— 


.04 
.03 





,03+° 


00 
.00 


WASP-3 


b 


1 


22 


185 -o.oo 


250 


9Q + 11.08 
-11.16 


0. 


07+.°, 


05 
.04 


, s , + ll.ll 
3-80_3 2i 


1 


72 


+0 
-0 


08 
.08 





,03« o 


00 

.00 


WASP-4 


b 





.91 


1 3 4+°-°° 
-0.00 


238 


gg + 10.83 
-11.83 


0. 


03+° 


04 
.02 


11 37+ 11 - 23 


1 


21 


+0 
-0 


06 
.06 





02+0 


00 

.00 


WASP-5 


b 


1 


.01 


1 63+°°° 
1 - bJ -o.oo 


282 


g2 +9.11 
3 -9.08 


0. 


07+° 


0-1 
.03 


9 77 +6.48 
^"-2.23 


1 


64 


+0 
-0 


.05 
.05 





,03+° 


00 
.00 


XO-1 


b 


1 


03 


3 94+°-°° 
J a -o.oo 


119 


n9 + 12.73 
U ^-11.60 


0. 


io+° 


11 

07 


3.09+ 8 / . 3 


0. 


94 


+0 
-0 


10 

.09 





0^° 


00 
.00 


XO-2 


b 





.97 


, ko+O.OO 

2 - b -o.oo 


86 - 00 ±io.07 


0. 


19+° 


15 
.13 


3 t.o + 9,'40 
• i - 68 -2.94 


0. 


56 


+0 

-0 


.06 
.07 







00 
.00 


XO-3 


b 


1 


41 


3 iq+0.00 

■ il9 -o.oo 




0. 


29+.°, 


01 

.01 


3Q 89 + 7 29 
Jaoa -6.39 


13 09+°- 14 
J -°- ua _0.14 





•°«ts 


00 
.00 


XO-4 


b 


1 


.32 


4 13+°-°° 
-0.00 


200 


, 9 +92.23 
-^-34.80 


0. 


35+.°, 


22 
.20 


4 34+13.O6 
-3.64 


1 


79 


+0 
-0 


50 
.27 





,06+0 


00 
00 
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Table 5. Table of the orbital parameters for a 2-planct fit, both directly output from EXOFIT and thence derived. The values of the 
parameters T, K, e and s (generated from exofit) are the medians of the parameter posterior distributions, with the associated 68.3% 
confidence regions. The other parameters were calculated using these values and stellar masses taken from the published literature. Note 
that some parameters are extremely well-constrained, hence the errors on the parameter estimates are so small as to appear to be zero 
to the two decimal places shown in this table. A full table in machine-readable format will be provided on the website, and the reader is 
directed there if such data arc required. 



planet 




m*(M soi ) 


T(days) 


Kims- 1 ) 




e 




s 






a(AU) 


BD-17 63 


b 


0.74 


1976 15+ 226 - 40 


297 6 2+ 389 - 34 
Z3,,tl -93 .47 





87+° 


07 
.07 


9 46+ 1 ' 96 
J ' 4 -1.47 


7 48+ 4 ' 67 
'■ 40 -1.37 


2 


.79+° 
-0 


21 
.17 










1 82 40+ 7 " 20 


0. 


84 +0 
o4 -0 


.01 
.01 




5 68+ ' 17 
°' DB _0.14 


3 


,54+ 9 
-0 


25 
.15 




b 


0.10 


1184 29 + 619 - 30 


870 93+ 6 00- 70 


0. 


43+° 
4J -0 


30 
.29 


29 49+ 159 46 

2J.42_ 27 44 


g i 5 +10.65 
a ' IS _6.11 


1 


02+° 
u -0 


33 
.97 








504 54+ 1270 37 
DU4 - D -493.52 


965 92+ 484 ' 80 
™ J ' 3 -952 .10 





39+° 


32 
.27 




6 84+ 9 ' 21 
°' M _6.80 





.58+° 
-0 


76 

.53 


silon Eri 
epsi on ri 


b 


0.86 


2443 29 + 45 ' 64 


13 qi+ 1 '^ 1 
U ' 91 -1.49 


0. 


05+° 


07 
.04 


_ 79 +0.84 
'''^-0.76 


83+ 09 
U ' S,3 -0.09 


3 


.37+° 
— 


04 
.04 








541 54+ 2 - 09 

D41.D4_ 1 Q5 


14 86+ 8 ' 97 
1 *- 0D _4.11 


0. 


82+° 


09 
.10 




30+ ' 08 


1 


.23+° 
••"'-0 


00 
.00 


c silon Tan 
cpsi on au 


b 


2.70 


bys - 14 _10 .56 


97 38+ 4 ' 08 
a '- JS -3.98 


0. 


13 +0 


06 
.05 


6 17+ 3 - 70 
C, ' i '-5.02 


7 76+ ' 30 
'•'°-0.31 


1 


93+° 
• a -0 


02 
.02 








71 16+1113.78 


12 48+ 13 ' 93 


0. 


66+° 


25 




18+ ' 90 





L, 

,47 + 


.58 








-70.55 


I ' i - 4S _10.77 


bt> -0 


.37 




"■ IS _0.14 


-0 


.45 




b 


1.59 


90^ 4^+ 4 - 49 


37 97+ 4 ' 81 
3'' y '_4.42 


0. 


04+° 
U4 -0 


07 

.03 


1 51+ 2 - 27 


2 46+ ' 31 
2i4b -0.29 


2 


.14+° 

-0 


.01 
.01 








Jb4s ' -661.35 


J_1 (in 

1761 05 + 17b - 92 
i(oi.uo_ 243 13 


0. 


75 + 


.09 




1 1 7 Q9 + 51 - 84 


5 


. . 4-2 

41 + „ 

-0 


08 
.68 


GJ3021 


b 


0.90 


133 71+ ' 20 


1 66 7^ + 3 96 


0. 


51+° 


02 
.02 


15 40+ 2 - 49 
10,4 -2.32 


, , K +0.08 
■ > - >O -0.08 





.49+° 


00 
.00 








oo (^o+3748.49 


5 00+ 14,75 
J ' UU -4.26 


0. 


56+° 


.38 




05 +0 ' 26 
Ui " J -0.05 





.20+ 4 


39 
.17 


GJ317 


b 


0.24 


4609 «q + 3780.22 
4bu ^ ba _ 1555.O8 


33 00 + 22 89 
*»-UU_ 7 .52 


0. 


35+° 
M -0 


.23 
.19 


4 51+ 4 ' 39 
4 ' bi -2.90 


95+ ' 86 
U ' ab -0.26 


3 


.37+ 1 
-0 


65 
.81 








682 64+ 4 ' 93 
DOZ - D -4.68 


82 68 + 4 ' 15 


0. 


27+° 


05 
.06 




, ,, + 0.06 

1 - > ' s -o.oe 





.94+° 


.00 
00 


GJ674 


b 


0.35 


34 82+ - 07 


4 84+°- 40 
4 - o4 _0.42 


0. 


23+° 
-"-0 


.08 
.08 


-1 , 7 +0.28 
1 ' 1 '-0.22 


04+ ' 00 
U ' U -0.00 





.15+° 

-a 


00 
.00 








4 69+ ' 00 


8 59+ 39 


0. 


14+0 


04 




03+ 00 





,04 +0 










4 ' ba -0.00 


SE>y -0.39 


i4 -0 


.05 






u -0 


.00 


GJ849 


b 


0.49 


1Q71 k-j + 178.78 


21 41+ 1 ' 79 
' !i - 4 -1.68 


0. 


11+° 


07 
.07 


2 89+ - 92 
^' SJ -0.85 


81+ ' 08 
U ' S1 -0.07 


2 


.43+° 
u -0 


.14 
09 








9065 cn+3863.21 


65 32+ 84 73 


0. 


59 +0 






3 06+ 2 95 


g 


.71+1 


79 








auoo.ou_ 3507 96 


DJ ' JZ -38.59 


OM -0 


.37 




a - UO -1.70 


— 1 


.87 


GJ86 


b 


0.80 


15 77 +o.oo 
lb -"-o.oo 


378 50+ 1-12 
■ 3 '°- ou -l.ll 


0. 


05+° 


00 
.00 


U ' MS -0.73 


4 02+ ' 01 
4,uz -o.oi 





.11+° 

-0 


.00 
.00 








6771 07+ 1223 " 25 
,1 - u ' -972.87 


1695 7l + 186 - 49 
id»o./i_ 272 g4 





80+° 


04 

.03 




80 82 + 7 " 66 


6 


50+° 


.76 
.64 


HAT-P-6 


b 


1.29 


87 07+ 4320 01 
8 '- u '_85.42 


„ „ 9 + 63.07 
8 - y2 -7.84 


0. 


55+° 


.33 
.38 


48 90 + 37 86 
48.2U_ 35 g9 


14+ 1 ' (i<i 
U ' 14 -0.12 





.42+ 5 
-0 


31 

.39 








, „, + 1126.88 

3 - 85 -o.oo 


1 07 1 1 + 1 3 - 37 
1U '- li _101.50 


0. 


11+° 


60 
.09 




1 00+ ' 14 
i ' UU -0.91 





.05+ 2 


.26 
.00 


HAT-P-8 


b 


1.28 


24 39 + 2 317-94 
- M - 3M _23.15 


,. 9 , + 49.73 


0. 


54+° 
°-0 


33 
.37 


4 94+ 4 - 87 
4 ' M4 -3.95 


13+ 1 ' 02 
"•"-O.ll 





.18+ 3 

-0 


57 
.15 








3 08+ ' 00 
3-US_ Q1 


160 51 + 809 

iou.ai_ 6 ?7 





05 +0 


06 
.03 




. ,_ + 0.07 
1 " 3o -0.06 





05+° 


00 
.00 


HAT-P-9 


b 


1.30 


, Q9 + 0.00 
3 a -0.00 


84 a<J +11 ' 66 


0. 


17+° 


14 

.10 


3 28+ 8 16 
J ' /8 -2.70 


77+ ' 10 
U '"-0.09 





.05+° 
— 


.00 
.00 








138 54 + 433 0- 13 
iJS - b4 -136 .58 


, q + 35.70 
8 - ia -7.18 


0. 


53+° 


33 
.37 




13+ 1 ' 31 





,57+ 5 
-0 


22 
.54 


HD101930 


b 


0.74 


362 8 2+ 3814 - 57 


2 85 + 7 S> ° 
^•°°-2.19 


0. 


55+° 


33 
.36 


1 62+ - 71 
-0.64 


05+ ' 26 
uuo _0.04 





90+ 3 
-0 


69 
.85 








70 64+ ' 48 


18 14+ ' 93 


0. 


08+° 


05 




30+ ' 02 





LO 

3<r 


00 








,U - M _0.44 


-0.92 


US -0 


.05 




U ' 3U _0.02 




u "-0 


.00 


HD108874 


b 


0.95 


39J.33_ Q g4 


37 76+ 1 14 
3''' t, _1.14 


0. 


07+° 


03 
.03 


81+ - 91 
U ' S1 -0.60 


1 ,,+0.04 
L - iL -0.04 


1 


.04+° 


00 
.00 








1 608 44+ 41 - 21 


18 4Q + 0-95 


0. 


OK + 






O qq+0 06 


2 


.64+° 


04 








1D0H.44_ 33 3g 


^■^-o.go 


J£l -o 


.04 




U ' aa -0.05 




u -0 


.04 


HD 11506 


b 


1.19 


1337 48+ 247 15 
1JJ '' 40 -64.76 


65 38+ 18 91 
b0JS -6.35 


0. 


30+° 


16 
.08 


4 08+ 1 ' 63 
4 ' US -1.35 


3 78+ 1 ' 15 
a ''°-0.35 


2 


,52+ 9 
-0 


.30 
.08 








1 70 '46 + 1 ' 82 
l(U.JD_ g9 93 


25 54+ 6 ' 16 





36+° 


13 
.17 




73+ - 16 
U ' '^-0.32 





.64+° 
u -0 


00 

.25 


HD118203 


b 


1.23 


6 13+ 000 

b -o.oo 


21 7 20 4 " 4 * 95 





31+° 
31 -0 


02 
.02 


14 ^4+ 3 - 29 

14.J4_ 2 g6 


2 14+ ' 04 
-0.04 





7+0 
' -0 


00 

.00 








6753 13+ 5243 ' 74 
DCJ3.13_ 3967 7£ . 


3J/.12_ 294 44 


0. 


51+° 


.24 

.33 




,, „ 9 +40.68 
31.82_ 24 Q2 


7 


,49+ 3 


50 
.34 


HD12661 


b 


1. 14 


969 67+0' 12 


74 Q-j + 0-73 


0. 


36+° 
■ }b -o 


01 
.01 


3 02+ - 67 
J ' U/ -0.64 


2 39+ ' 02 





.84+° 
-0 


00 
.00 








1681 47+ 29 - 15 


29 16+ ' 79 


0. 


02+° 


02 




1 86+ 05 


2 


,89+ 9 


03 








lbB1 ' 4 '-26.29 


^ M ' 1D _0.83 


U ^-0 


.01 




1 ' SO -0.05 






HD128311 


b 


0.83 


924 69+ 6 43 
BZ, -° 3 -6 .83 


75 4S + 3 ' 28 


0. 


14+0 


10 
.10 


16 1 c^ + l-- 55 
- LO -- LJ -1.37 


3 14+0- 12 
a ' 1 ^-0.12 


1 


■ 74 1°0 


01 

.01 




c 




458.32+ 3 09 
— 2.95 


53.71+=- 39 


0. 


33+0 
— 


06 
.06 




1.69+ 9 '! 8 
— 0.19 


1 


.09+. 9 


.00 
.00 


HD131664 


b 


1.10 


1964 5i+ 35 - 83 
ab b -45.84 


363 85+ 31 98 


0. 


b4 l°o 


03 
.03 


4 67+ - 86 
-0.74 


IS 9S+1-10 
18 - 28 _0.86 


3 


"IS 


04 

.05 




c 




c«i ,., + 1790.01 
561 - 36 -558.23 


5 .4o; 4 ° 2 3 9 2 


0. 


«l°o 


39 
.31 




14+ ' 53 
° -0.13 


1 


,37+ 2 


.20 
.33 


HD132406 


b 


1.09 


540 9l+ 892 18 
ow.ai_ 32g 23 


49.331, 8 ; 99 


0. 


.391° 


.31 
.24 


6.90+ 9 73 


1 74+4-10 
-1.19 


1 


,34+J 


.22 
.62 




c 








0. 


.391° 


.29 
.25 




4 75+ 2 ' 15 
-3.97 


2 


.22+° 


46 

.23 


HD142 


b 


1.23 


349.57+. 3 ; 49 


29.66; 4 ;f 2 


0. 


20+J 


11 
.11 


8.74+i;»i 


1 15+ - 16 
ii5 -0.15 


1 


° 4 1°0 


.01 
.01 




c 




9899 QS + 3348 - 44 
9S22.95_ 320g 17 


t>3.42_ 17 5g 


0. 


18l°0 


.20 
.13 




6.27; 7 ;J 3 


9 


,62+ 2 


08 
.23 


HD142022 


b 


0.90 


244 04+ 6905 ' 04 
u -242.28 


3 41+ 861 
-2.65 


0. 


™ + -°0 


.35 
.36 


2.84+ i; ™ 


05+ - 52 
U ' UE> -0.05 





74 l 6 o 


27 
.71 




c 




1877.47+ 3 ^; 74 


127 61+ 84 60 


0. 


62+ 9 


11 

.10 




5 70+ 2 - 52 
°-' U -1.21 


2 


.8 8 ; 9 


.04 
02 


HD149143 


b 


1.20 


4 07 +0.00 

-0.00 


149.85 + V 79 
— 1 .66 


0. 


0ll°0 


01 
.00 


, ,, + 1-83 


, ,, + 0.02 
i - 3J -0.01 





■"IS 


00 
.00 




c 




248 47+ 33 °0.06 


, 7 , + 13.28 


0. 




33 




07+ ' 55 





K9 + 4 

■ 82 -o 


02 








™ ,4 '-245.40 


-3.00 




.37 




-0.06 




.78 


HD154345 


b 


0.89 


3216.13; 3 f 70 8 | 3 




0. 


10l°o 


47 

.08 


9 ,„+0.47 
2 ' 38 -0.41 


„ S7 +0.12 
°' 87 -0.84 


4 


.10+ 9 


13 
.86 




c 




2695.22t 2 97 g 19 6 


6.21+ 8 ; 39 


0. 


15l°0 


55 
.11 




15+ - 83 
U ' 1J -0.13 


3 


,65+ 9 


.60 
.41 


HD155358 


b 


0.87 


193. 24+ 9 ; 98 


32.88+J;^ 


0. 




04 
.04 




O 84+ - 04 
°- 84 _0.04 





.62+_° 


.00 
00 




c 




314.82+ 4 ; 92 


33 fi -) + 41 -74 
33.03_ 12 ?7 


0. 


83+° 


12 
.17 




,, ,7+0.23 
U ' j7 -0.11 





.86±° 


01 

.01 


HD162020 


b 


0.80 


o ,,+0.00 

8 - 43 -o.oo 


1808.83+ 5 5 | 


0. 


28i° 


00 
.00 


io.«+«; 


15 00 + 04 
15 ' UU _0.04 





.08±° 


.00 
.00 




c 




180 71 + 32 7159 
180.71_ 1?5 94 


6.73+JV/ 


0. 


S 7 1°0 


.33 
.39 




io+°' 41 

u ' lu -0.09 





™tl 


57 
.53 


HD168443 


b 


1.01 


58 11+0-00 


529.04+ 82 ; 29 


0. 


"l°o 


04 
.04 


7.80t 3 ; 74 


„ , 7 +0.93 
8 ' 37 -0.61 





.29+_° 


00 
.00 




c 




1755.58+.^ 


302. 53+ 3 ' 8 ^ 


0. 


23+_° 


02 
.02 




l761 +0.21 
-0.21 


2 


.86t° 


.01 
.01 
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A Uniformly Derived Catalogue of Exoplanets from Radial Velocities 



planet 



*.(M so! ) 



Kims- 1 ) 



a(AU) 



HD169830 b 

c 

HD171028 b 

c 

HD183263 b 
c 

HD185269 b 
c 

HD187123 b 

c 

HD189733 b 

c 

HD190228 b 

c 

HD190360 b 

c 

HD190647 b 

c 

HD195019 b 
c 

HD202206 b 

c 

HD20868 b 
c 

HD209458 b 

c 

HD212301 b 

c 

HD217107 b 

c 

HD219828 b 

c 

HD221287 b 
c 

HD224693 b 
c 

HD23127 b 
c 

HD2638 b 
c 

HD27442 b 
c 

HD27894 b 
c 

HD28185 b 
c 

HD285968 b 
c 

HD330075 b 
c 

HD33636 b 
c 

HD3651 b 

c 

HD38529 b 
c 

HD4203 b 
c 

HD4208 b 
c 

HD43691 b 
c 

HD43848 b 
c 

HD46375 b 



1.12 
1.30 
1.04 
0.81 
1.82 
0.98 
1.10 
1.02 
1.07 
0.78 
1.13 
1.05 
1.11 
1.24 
1.30 
1.30 
1.13 
0.93 
1.20 
0.75 
0.99 
0.49 
0.70 
1.02 
0.88 
1.48 
1.13 
0.88 
1.38 
0.93 
0.93 



(i.84 1 



ook ri +0.29 
225 - 61 -0.31 
13Q eo + 4188.83 

498 - 38 ±4i:o2 

267.54+. 7 ; 3 2 

4830.04+i«| : ™ 
A +0.00 
-0.00 
o, Q „ + 2639.27 
32a8 _31.94 

^toTo 
5502.18+.™ 
2 22 +0.00 



-0.00 
149 82 + 3434 - 51 
1 * a - ! "-147 .94 

1140 04+ 16 47 
-15.62 

2925.88+ 3, J 4 | 

184.70+™ 

1036.28+» 4 27 , 5 

422 54+3!3 4 .9 4 

45 - 87 _0.76 
51 93+ 08 
7 +0.08 



255. 87 1 
258.20 



0.08 
+ 1.14 



380.862 

m.osj 

3.52+ 



-1.02 
, + 0.08 
-0.08 
- + 1942.40 
-110.42 
,+0.00 
-0.00 
OR/I fi , + 1050. 15 
264 - 63 -259.91 

R1 +6274.66 
2356.61_ 21Q4 ?0 

, ,,.+0.00 
^ D -0.00 
- , ,+0.00 
'■ i3 -0.00 

4106.23+- 8 3.° 8 5 8 

3.83«» 

956.74;--/ 

°- 60 -o.oo 

455 12 + 6 ' 45 

96 7o+°- 03 
2b - 7o -0.03 

7882.48; 42 f 9 5 2 5 

1124.66^t^i 10 44 
1219 i4+ 36 - 54 

-1099.09 

3 44+° 00 
-0.00 

26 04+ 1978 12 
^ u ' u -25.00 

410 58 + 866 97 
i±u.oo_ 407 42 

417 14+9 4 5-92 

23.31+. 9 - 

17.98+» ° 3 

381.37+1; 22 

508 nq+ 1160 - 41 
508. 09_ lg8 94 

10.2 5 ; 8 4 | 8 7 - 88 

10 23+ 01 

, „„+0.00 



655.66 



+ 6970-11 
-589.87 
-,+11.43 



2 l 27 -7lll . 98 
7K41 o, + 3695. 12 
7841.32_ 3270 74 

294 K7+ f > 1 - 82 
3 "'-195.79 

62 25+ 03 

1431 +o.oo 

° -0.00 

2148 41+ 595 
-7.87 

« 8 -0 4 iI'69 

O + 202.96 
-200.21 

129 04+ 708 - 37 
au -121.12 

829.27^ 9 -= 7 

36.96+.° ° 4 

22 6 1 - 82 1 5 22 3 52 9 19 
2305.82+f 2 8 61 97 4 
o,+0.00 



391.! 



83.01+. 3 ; 29 

3.39t»;« 

72.24±%% 
47.44±«;« 
86.82±\- 2 2 l 
74 41 + 4 2.59 
-19.98 
89.8ll 4 ; 3 | 

4 -"i 9 3 :f6 

28.50+J; 99 
R +4.00 
-4.07 



198. 36"* 



25.58 
2.23 



+ 4.20 
■3.99 
+ 5.57 



-1.79 
92.35^; 49 

n +10.61 
-4.14 
7+189 



5.10"' 



159.1 



36 ' 17 -1.75 
2.88l 2 3 2 f 
17R S7+70.77 
176 - 87 -43.26 
a + 33.08 
-34.85 

706.93 + " 9 „ 27 

— 65.64 
291 17+16 4 -81 

100.22 + ° ^| 

1 59+ 219 

84.27+° ; 8 4 
9 s r + 3.04 
2 - 85 -2.03 



26.81 



+ 76.37 



56. 85^ 



-17.28 
- + 1.92 
-1.93 
138.34+111 
39 48+ 13 * 92 
7 99 + 0.59 
'• -0.62 
78 15+ 13 !- 30 
' SilJ -45.22 

9.66+ 6 8 ; 3 i 

71.62+ 7 6 1 f 

38.90+J 77 

1849 37+ 94 - 8 l 
iMa -' 3 '-501.16 
o ,o + 19.55 
8 - 38 -7.33 

4 - 4 ^ 9 2 : 2 3°6 

O + 26.90 
-5.26 
g+2.69 
-28.78 
10.37+ 2 - 8! J 

56 34+ 2 -56 
Jt> ' 3 -1.29 
161 1,+ 8 -85 

22 84+ 26 ' 97 
"-11.21 

3 - 84 i3.20 
10.62+ 2 . , 9 

106.96+1 ™ 

10 ,i+ 98 -' 32 
iU - 0i -8.53 
770 s . + 123.39 
773.84_ 7f . 2 03 

1923. 73+?!, 92 



6.33] 
31.26"" 



Q +l.i 



-174.57 



3.02; 
30.41 



0.00 
+ 931.01 



3 ' 49 ll.60 

54.97+-™ 

170.83+1; 92 

56.24+ 29 5 9 l 

10 70+11- 08 

2 62+ 3 03 
^■'"-l.so 

18.791 9 ; 82 
19.26+.™; 84 
125.00+. 4 ; 49 
44 31+612-06 
-41-42 

476 62+ 555 76 
«' D - oz _471.61 

33 79+ - 79 
3J ' a -0.81 



„ O-+0.03 

U ' 3 '-0.03 

55 + - 33 
U55 -0.37 

61+ ' 16 
° -0.11 

57+°-l 3 
U ' a7 -0.08 

n , 0+0.01 

°' 38 -0.01 

07+ - 09 
u -0.05 
n ,o + 0.03 
°- 28 _0.04 
55+ - 32 
°' 55 -0.36 

01 + °1 
-0.01 

32+ 12 

02 +0 02 

96+° °l 
-0.01 

55+ - 33 
U ' a5 -0.37 

0.52+°; 9 4 
0.34+°;" 

55+ - 32 
°' 55 -0.35 

17+ 003 
" ' -0.04 
„ ^, + 0.30 
°' 56 -0.36 
a+0.11 



79" 1 
u ' ' M -0.17 

74+ - 07 
-0.08 
„ O-+0.02 
U " rl '-0.02 

26 + 06 
U ' 2b -0.06 

76+ - 00 
-0.00 
n +0.30 
-0.38 

01 + °1 

" ° -0.01 

54+ - 34 
u -0.37 
,+0.35 
0.36 
+ 0.04 
-0.04 
+ 0.01 
-0.01 
0.04 
0.03 
+ 0.09 
-0.06 
0.17 
0.20 
+ 0.39 
0.28 
+ 0.12 
0.08 

o-o 4 i°o:°o 4 3 
0.68+s- 

46+ - 34 
" "-0.24 

44+°-l 7 
" -0.14 

01 + °1 
-0.01 
-, +0.39 



0.60"' 



0.52 
0.09 
0.13 
0.56 
0.09 
0.53 
0.41 
0.13 



0.41" 



2.93 



+ 4.85 
-2.28 



0.27 
-, + 0.52 



1()T"-"" 
-0.06 
„ ,o + 0.44 
°' 28 -0.22 

05 + 03 
U ' Ut> -0.03 

04+ 03 
-0.03 
„ 37 +0.36 
U ' 3 '-0.21 

43+°- 4 l 
-0.32 
q+0.29 
0.14 
+ 0.01 
-0.00 
+ 0.40 
-0.33 
+ 0.03 
-0.22 
+ 0-03 
0.01 
+ 0.41 
-0.23 

60+° 05 
u ' DU -0.06 

17+ 003 
-0.03 

34+°-°l 
-0.01 

69+° l 3 
U ' Oa -0.08 
,+0.33 
-0.23 
53+ ' 34 

06+ - 05 
-0.04 

45+ - 36 
-0.31 

12+ ' 04 
U ' i2 -0.04 
7+0.16 
0.37 
73+° l 2 
U73 -0.35 
„ nc - + 0.03 
°' U5 -0.03 
e; + 0.25 
0.42 



0.19 
0.01 
0.42 
0. 

0.61 
0.32 



0.32"' 



0.67 1 



0.65 1 



-1+1-64 
-2.14 



149.9 ; 9 ; 7 7 



74+ - 43 
-0.48 



io.i9t 9 ; 6 59 



, 70+0.55 

l- 78 -0.39 



11 09 + 3 ' 45 
ll' Ua -2.98 



2.18 
+0.04 



2 91+°-ll 
" -0.11 
„ nc -+0-43 
U - U5 -0.05 
2 27 +1.75 
-0.69 
-, ,„+0.18 

o RC -+0-05 
3 ' b5 -0.05 
6.S7+ 472 

°' 96 -0 + 

05 +0 -l 9 
-0.04 
„ e--, +0.00 
u -0.00 

1-"1 0'.°0 2 2 

05+°-°l 
-0.01 

05+ - 30 

6 07+ 4 92 
" -0.64 
-, ,o+0.16 

09 +0 49 
U ' Ua -0.08 

1 89+°-l° 
1 ' Ba _0.10 
„ „o+0.77 
°- 08 -0.07 

1 96+ - 64 
aD -0.42 

1 95+ - 52 
!' aD -0.45 

21 46+ 3 - 75 
"-2.07 

„ ,0+4.97 

y ' 23 -2.03 

1 99+ °l 
1 - aa _o.oi 

0.02+« « 7 

o.68«;«i 

07+° ■ 11 
-0.06 
-, ,,+4.93 
!^^-0.89 
„ oo+O.Ol 

O-^S-o.oi 

3+0.01 



i-SSlo.Ol 
„ 77+1.00 
^■"-0.27 

06+°-°l 
-0.01 
r + 10.00 



3 - 55 l 2 .29 
05 +0 -ll 

,-, +0.65 
3 ' 2i -0.32 
„ K o+0.03 
0<38 -0.03 

1 53 4S + 42 - 92 
15 3- 4S _58.00 

H+ L28 
°' 15 _0.14 

1 3-1+O.23 
l' 3i -1.26 
n 4K+° °1 
°' 48 -0.01 

t;+0.34 
0.03 
+ 1.23 
0.12 
1.29+0-11 

o.i2+" ;o 2 

o.60«;« 3 

-, +0.32 
0.38 

o- 8 ^ :°49 

03+ 05 
UU3 -0.02 
7+0.01 



0.05" 
0.14 



5.71" 



°- 07 Io,04 
62+°-°l 

0.27+ 5 27 
22 - 86 t22 3 32 



149.03" 



0.10 
0.23" 



+0.05 



27.71 



0.02 

°' 84 -0.02 

13 23+°-l 4 
i3 - /3 -0.14 
-, 7 ,+0-53 
-0.38 

37+ - 42 
UJ '-0.20 

o-o^° ;° I 

o-^rot 

„ e-7+5.63 
°- j7 -0.39 
, ,,+0.10 
2 -' ,2 -0.10 
j 41 + 26.25 

22 i6+i 3S 50 
^"•- Lu -22.11 

23+°-°l 
U - /3 -0.01 
„ Qo+0.08 
U - U,i -0.02 

. . . continued < 



n 01 +0.00 
°' 81 -0.00 
59+ 5 - 24 
°' 5a -0.56 
O + 0.05 
-0.07 
+ 0.01 
0.01 
+ 0.00 
-0.00 



1.23"* 
0.81 



5 K^ti-w 

5 - t,5 _0.91 
„ oo + O.OO 
°- 08 -0.00 

U '^ 2 -0.20 

04+ 00 
-0.00 
R ,o + l. 96 

6 - 18 -1.03 
+ 0.00 

-0.00 



0.03"* 



05+ 00 
U ' Ut> -0.00 

67+ 4 - 92 
u -0.64 

2 61+ - 03 
2 - 6i -0.02 

qo + 0.03 

+ 4.30 
-0.59 

2 07+°-°l 
ZM -0.02 

1 14+ 3 -57 
-1.04 

- + 0.04 
-0.00 
7+0.00 
-0.00 
-, +0.00 
-0.00 



0.63^ 



0.25 

o.2r<_ 

0.81" 1 



0.95" 1 



0.03 1 



1.26 



n 01+O.OO 
°- 81 -o.oo 

^+0.00 

-0.00 

42+ 2 - 49 
-0.40 

05 + 00 

u - ua -o.oo 

3.52+ 4 ; 8 / 3 
,+0.00 
-0.00 

„ 00+0.00 
°- 08 -o.oo 

5 20+°- 2 l 
5 - 20 -0.10 

05+ 00 
u - u& -o.oo 

2 04+1" 58 
^■ U *-0.56 

02+1- 38 
UU ^-0.00 
+0.01 
0.01 

0- ^toTo 
8.46+ 2 ;f 5 

, ,„+0.99 
Ji!U -2.16 

2.33+°; 9 / 3 

04+ - 00 

-0.00 
12.87 
0.15 
t;+1.30 
-1.10 

1 16+1- 40 

-0.02 

15+ - 04 
U 'l b -0.00 

12+ - 00 
U ' i2 -0.00 

,+0.00 
0.00 

1 24+1- 50 

-0.35 

07+1- 92 
-0.05 

07+ - 00 
u -0.01 

4+0.00 
-0.00 

1 31 + 5 - 42 

1- 3i -1.03 

3 2K+0-01 

7.78+ 2 ; 28 
o.83+S;ii 

„ ,„+0.00 

°- 2a -o.oo 

13+°-°° 

u - i3 -0.00 

3.71" 1 



0.17] 
1.15+ 



1.03" 



0.04"f 



1.18 
1.09 



+0.01 
0.01 
+ 0.01 
0.01 
+ 0.35 
0.41 
5+1.19 



°- 48 lo,40 

1 66+ °l 
-0.01 

i.40+?; 74 3 
O-^oTo 

3.29+0- 48 

o O0 + 0.15 

d - di -3.10 

04+ 00 
-0.00 

19+1- 67 
u - la _0.16 
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planet 




TTT *. ( M i 1 
'"■*\ 1VI SOI ) 


T(days) 


K(ma — 1 ) 










m p (M Jup ) 




a(AU) 


HD47536 


b 


0.94 


698. 09+ 24 ,, 29 

— 686. !JG 


10512 t 2 98.90 





18 i° 


14 


4 58+ 9 " 07 
J -3.83 


4 


m Q5 

3 3 ■ ^ 
°°-4.28 


1 


51 i 1.41 




c 




713 - 17 i 6 62 35 


87.58+4 1 i ; 2 4 





solo 1 


20 




3 


«^+l 67 
76 

-3.69 


1 


53 -1.27 


HD49674 


b 


1.01 


275.09+.^° 


4.23+. 4 2 j 7 







42 


2.56+!" 94 





. . + 10 

10^ 

-0.07 





oo + 3 .33 
83 -0.65 




c 




4 - 95 loTo 


11.8811;"; 





09 i° 


06 







. ^+0 01 
10^ 

-0.01 





06 +o.oo 

u -0.00 


HD50499 


b 


1.28 


2450.67+., 4 ; 98 


21.82+1;^ 





2 si° 


07 


2 ' 95 -1.08 


1 


,. 4-0 1 5 
b5 -0.15 


3 


of. + 0.06 
86 -0.05 




c 




9734.20t| 4 66.65 


57 00+ 109 " 73 
°'' UU -32.17 





5 4 lo 


34 




5 


b5 -2.96 


9 


gq + 2.17 
° -2.46 


HD5319 


b 


1.60 


684.85+^.55 


37.35+ 9 ' 7 2 





iiio 




5- 43 to.82 


2 


zu -0.28 


1 


V o + .02 
78 -0.04 




c 




1541.64+^f 


15.24+™» 





4 ll°o 


29 
28 




1 


^,-,-t-O 56 

08+0;^ 


3 


ne + 1.93 
U1> -0.57 


HD63454 


b 


0.80 


j_ qoor; rii 

iiLL -32.19 


3 - 07 ±2.50 







39 


5.54+1; 3 ! 





,,^. + 15 
03 ±0.02 





M -0.17 




c 




2.82+« ; oo 


63.26+.1; 78 





02+ 9 









rt ,-. + 01 

38 Io.oi 





04+ 00 














02 






u -0.00 


HD689S8 


b 


1.12 




189.5711^ 







01 


3.36^! 1? 


1 


-fO.Ol 
83 -0.01 





q7+0 ■ 00 
— 0.00 




c 






68.13™ 





16lo 






5 


^«4-l 68 
° y -0-72 


5 


17+1.14 
-0.51 


HD73267 


b 


0.89 


1265 ' 34 ±9.40 


18.63™ 





56i°0 


2 7 

53 


-"tS:iS 





724 : .43 


2 


,n+0-01 
^ U -0.01 




c 




1278.04+.1S; 28 


73.21^ 9 3 9 3 7 





28 i° 


09 




3 


.^4-1 00 
^"-0.23 


2 


fiii^O . 02 
-0.03 


HD73526 


b 


1.01 




45 - 30+ .6'41 





7 °l°o 


09 


9 - 5 ^2.23 





+ 0.13 
91 -0.11 





(^ry-\-0 . 00 

b -0.00 




c 




193.42+J;« 


114 30+ 8 ' 71 
-8.35 





5 2 i° 






2 


7( i+0.16 
-0.16 





t;^ + .00 
66 -0.00 


HD74156 


b 


1.24 


B1 - aB ±S:Si 
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Table 6. The number of published planets compared with the best-fit model from this analysis (i.e. the flags and reduced chi-square ratios 
for the results for each system). The 'candidates' column shows the current number of confirmed planets (from http://www.exoplanet.eu 
and http://exoplanets.org, as of 201f August 01), and the 'visual quality flag' (assigned by eye) is the best EXOFIT model, where '1' 
signifies that the 1-planet fit is best, and '2' means that the 2-planet fit is best. '3' means that both 1- and 2-planet solutions provide 
equally good or bad fits, and this class is again subdivided into '3a' and '3b', as explained in Section 6.2. Also shown is the log-likelihood 
ratio of the chi-square values, R, as defined in Section 5. Our visual flag assignments are validated somewhat by noting that in 99% of 
systems the visual flag and chi-square results agree (or at least are not contradictory, for the class 3 cases). Those systems denoted by '-' 
are those where there were not sufficient degrees of freedom to calculate a value for the log likelihood ratio. The prior flag is also shown, 
where flag TV indicates that the analysis was performed using the normal priors shown in Tables 1 and 2, and flag D indicates an analysis 
with different priors as shown in Table 3. 



system number of candidates 
from literature 


visual quality flag 


R 


flag for period 
prior used 


BD-17 63 


1 


1 


25.06 


N 


ChaHa8 


1 


1 


- 


N 


cpsilon Eri 


1 


3a 


-0.84 


D 


cpsilon Tau 


1 


1 


4.62 


N 


gamma Ccp 


1 


2 


-398.83 


D 


GJ3021 


1 


3b 


0.28 


N 


GJ317 


1 


1 


47.75 


N 


GJ674 


1 


3a 


-8.44 


N 


GJ849 


1 


3a 


8.20 


D 


GJ86 


1 


2 


-1554.76 


D 


HAT-P-6 


1 


1 


1450.29 


N 


HAT-P-8 


1 


3b 


- 


N 


HAT-P-9 


1 


3b 


0.95 


D 


HD101930 


1 


3b 


6.12 


N 


HD108874 


2 


2 


-8.23 


N 


HD11506 


2 


2 


-3.36 


N 


HD118203 


1 


1 


80.12 


D 


HD12661 


2 


2 


-15.22 


D 


HD128311 


2 


2 


-32.30 


D 


HD131664 


1 


1 


3.34 


D 


HD132406 


1 


1 


234.41 


N 


HD142 


1 


3a 


-5.87 


D 


HD142022 


1 


3b 


0.43 


N 


HD149143 


1 


3b 


1.49 


N 


HD154345 


1 


1 


44.33 


N 


HD155358 


2 


2 


-2.90 


N 


HD162020 


1 


3b 


102.59 


D 


HD168443 


2 


2 


-994.46 


D 


HD169830 


2 


3b 


0.25 


D 


HD171028 


1 


3a 


18.72 


N 


HD183263 


2 


2 


-187.35 


D 


HD185269 


1 


3b 


0.83 


N 


HD187123 


2 


3b 


-29.06 


D 


HD189733 


1 


3a 


-140.98 


D 


HD190228 


1 


3b 


0.08 


N 


HD190360 


2 


3b 


0.04 


D 


HD190647 


1 


3b 


5.01 


N 


HD195019 


1 


1 


771.62 


N 


HD202206 


2 


3b 


-4.59 


D 


HD20868 


1 


2 


-252.94 


D 


HD209458 


1 


3b 


-0.03 


D 


HD212301 


1 


1 


5.15 


N 


HD217107 


2 


3a 


-59.88 


D 


HD219828 


1 


3a 


-352.73 


D 


HD221287 


1 


1 


141.08 


N 


HD224693 


1 


3b 


20.92 


N 


HD23127 


1 


1 


69.47 


N 


HD2638 


1 


3b 


-0.97 


N 


HD27442 


1 


3b 


10.37 


N 


HD27894 


1 


1 


53.60 


N 


HD28185 


1 


1 


10.27 


D 


HD285968 


1 


1 


0.42 


N 


HD330075 


1 


3b 


745.43 


D 


HD33636 


1 


1 


104.05 


D 


HD3651 


1 


3b 


-0.06 


N 


HD38529 


2 


2 


-41.83 


D 


HD4203 


1 


3a 


4.92 


N 


HD4208 


1 


3b 


0.65 


N 


HD43691 


1 


3b 


5.07 


N 


HD43848 


1 


1 




N 


HD46375 


1 


3b 


-0.04 


D 


HD47536 


2 


1 


17.15 


D 


HD49674 


1 


3b 


0.24 


N 


HD50499 


1 


2 


-6.74 


D 
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system 


number of candidates 
from literature 


visual quality flag 


R 


flag for period 
prior used 


HD5319 


1 


2 


-12.72 


Y) 


HD63454 


1 


3b 


2.57 


N 


HD68988 


1 


2 


-100.51 


Y) 


HD73267 


1 


3b 


7.60 




HD73526 


2 


2 


-7.43 


N 


HD74156 


2 


2 


-241.14 


-£) 




I 


I 


0.14 


]\ 




1 


3a 


28.79 


]\* 


HD76700 


I 


1 


14.28 


N 


HD80606 


1 


I 


4.22 


-£) 


HD81040 


I 


3b 


-1.13 




HD82943 


2 


2 


-319.74 


N 


HD8574 


1 


I 


-0.03 


]\- 


HD86081 


1 


x 


0.19 


-£) 


HD89307 


1 


3b 




-£) 


HD93083 


1 


3b 


32.00 




HR 810 


1 


3b 


8.71 


N 


kappa CrB 


1 


1 


8.87 


N 


NGC 2423 3 


1 


3b 


52.75 


N 


NGC4 349 127 


1 


3b 


24.85 


N 


tau Boo 


1 


3b 


-8.90 


D 


TrES-3 


1 


2 




D 


WASP-2 


1 


3b 




D 


WASP-3 


1 


3b 




D 


WASP-4 


1 


3b 


-33.89 


N 


WASP-5 


1 


3b 




N 


XO-1 


1 


1 




D 


XO-2 


1 


3b 




D 


XO-3 


1 


1 


20.24 


N 


XO-4 


1 


3b 




D 



